Nanopore Guided Annotation of Transcriptome Architectures.

High-resolution annotations of transcriptomes from all domains of life are essential for many sequencing-based RNA analyses, including Nanopore direct RNA sequencing (DRS), which would otherwise be hindered by misalignments and other analysis artefacts. DRS allows the capture and full-length sequencing of native RNAs, without recoding or amplification bias, and resulting data may be interrogated to define the identity and location of chemically modified ribonucleotides, as well as the length of poly(A) tails on individual RNA molecules. Existing software solutions for generating high-resolution transcriptome annotations are poorly suited to small gene dense organisms such as viruses due to the challenge of identifying distinct transcript isoforms where alternative splicing and overlapping RNAs are prevalent. To resolve this, we identified key characteristics of DRS datasets and developed a novel approach to transcriptome. We demonstrate, using a combination of synthetic and original datasets, that our novel approach yields a high level of precision and recall when reconstructing both gene sparse and gene dense transcriptomes from DRS datasets. We further apply this approach to generate a new high resolution transcriptome annotation of the neglected pathogen human adenovirus type F 41 for which we identify 77 distinct transcripts encoding at least 23 different proteins.

Novel viral splicing events and open reading frames revealed by long-read direct RNA sequencing of adenovirus transcripts.

Adenovirus is a common human pathogen that relies on host cell processes for transcription and processing of viral RNA and protein production. Although adenoviral promoters, splice junctions, and polyadenylation sites have been characterized using low-throughput biochemical techniques or short read cDNA-based sequencing, these technologies do not fully capture the complexity of the adenoviral transcriptome. By combining Illumina short-read and nanopore long-read direct RNA sequencing approaches, we mapped transcription start sites and RNA cleavage and polyadenylation sites across the adenovirus genome. In addition to confirming the known canonical viral early and late RNA cassettes, our analysis of splice junctions within long RNA reads revealed an additional 35 novel viral transcripts that meet stringent criteria for expression. These RNAs include fourteen new splice junctions which lead to expression of canonical open reading frames (ORFs), six novel ORF-containing transcripts, and 15 transcripts encoding for messages that could alter protein functions through truncation or fusion of canonical ORFs. In addition, we detect RNAs that bypass canonical cleavage sites and generate potential chimeric proteins by linking distinct gene transcription units. Among these chimeric proteins we detected an evolutionarily conserved protein containing the N-terminus of E4orf6 fused to the downstream DBP/E2A ORF. Loss of this novel protein, E4orf6/DBP, was associated with aberrant viral replication center morphology and poor viral spread. Our work highlights how long-read sequencing technologies combined with mass spectrometry can reveal further complexity within viral transcriptomes and resulting proteomes.

Epstein-Barr Virus-Encoded BILF1 Orthologues From Porcine Lymphotropic Herpesviruses Display Common Molecular Functionality.

Infection of immunosuppressed transplant patients with the human γ-herpesvirus Epstein-Barr virus (EBV) is associated with post-transplant lymphoproliferative disease (PTLD), an often fatal complication. Immunosuppressed miniature pigs infected with γ-herpesvirus porcine lymphotropic herpesvirus 1 (PLHV1) develop a similar disease, identifying pigs as a potential preclinical model for PTLD in humans. BILF1 is a G protein-coupled receptor (GPCR) encoded by EBV with constitutive activity linked to tumorigenesis and immunoevasive function downregulating MHC-I. In the present study, we compared BILF1-orthologues encoded by the three known PLHVs (PLHV1-3) with EBV-BILF1 to determine pharmacological suitability of BILF1 orthologues as model system to study EBV-BILF1 druggability. Cell surface localization, constitutive internalization, and MHC-I downregulation as well as membrane proximal constitutive Gαi signaling patterns were conserved across all BILFs. Only subtle differences between the individual BILFs were observed in downstream transcription factor activation. Using Illumina sequencing, PLHV1 was observed in lymphatic tissue from PTLD-diseased, but not non-diseased pigs. Importantly, these tissues showed enhanced expression of PLHV1-BILF1 supporting its involvement in PTLD infection.

Single-cell transcriptomics identifies Gadd45b as a regulator of herpesvirus-reactivating neurons.

Single-cell RNA sequencing (scRNA-seq) is a powerful technique for dissecting the complexity of normal and diseased tissues, enabling characterization of cell diversity and heterogeneous phenotypic states in unprecedented detail. However, this technology has been underutilized for exploring the interactions between the host cell and viral pathogens in latently infected cells. Herein, we use scRNA-seq and single-molecule sensitivity fluorescent in situ hybridization (smFISH) technologies to investigate host single-cell transcriptome changes upon the reactivation of a human neurotropic virus, herpes simplex virus-1 (HSV-1). We identify the stress sensor growth arrest and DNA damage-inducible 45 beta (Gadd45b) as a critical antiviral host factor that regulates HSV-1 reactivation events in a subpopulation of latently infected primary neurons. We show that distinct subcellular localization of Gadd45b correlates with the viral late gene expression program, as well as the expression of the viral transcription factor, ICP4. We propose that a hallmark of a "successful" or "aborted" HSV-1 reactivation state in primary neurons is determined by a unique subcellular localization signature of the stress sensor Gadd45b.

Targeting the m6A RNA modification pathway blocks SARS-CoV-2 and HCoV-OC43 replication.

N6-methyladenosine (m6A) is an abundant internal RNA modification, influencing transcript fate and function in uninfected and virus-infected cells. Installation of m6A by the nuclear RNA methyltransferase METTL3 occurs cotranscriptionally; however, the genomes of some cytoplasmic RNA viruses are also m6A-modified. How the cellular m6A modification machinery impacts coronavirus replication, which occurs exclusively in the cytoplasm, is unknown. Here we show that replication of SARS-CoV-2, the agent responsible for the COVID-19 pandemic, and a seasonal human ß-coronavirus HCoV-OC43, can be suppressed by depletion of METTL3 or cytoplasmic m6A reader proteins YTHDF1 and YTHDF3 and by a highly specific small molecule METTL3 inhibitor. Reduction of infectious titer correlates with decreased synthesis of viral RNAs and the essential nucleocapsid (N) protein. Sites of m6A modification on genomic and subgenomic RNAs of both viruses were mapped by methylated RNA immunoprecipitation sequencing (meRIP-seq). Levels of host factors involved in m6A installation, removal, and recognition were unchanged by HCoV-OC43 infection; however, nuclear localization of METTL3 and cytoplasmic m6A readers YTHDF1 and YTHDF2 increased. This establishes that coronavirus RNAs are m6A-modified and host m6A pathway components control ß-coronavirus replication. Moreover, it illustrates the therapeutic potential of targeting the m6A pathway to restrict coronavirus reproduction.

The influence of human genetic variation on Epstein-Barr virus sequence diversity.

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm3 and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

Direct RNA sequencing reveals m6A modifications on adenovirus RNA are necessary for efficient splicing.

Adenovirus is a nuclear replicating DNA virus reliant on host RNA processing machinery. Processing and metabolism of cellular RNAs can be regulated by METTL3, which catalyzes the addition of N6-methyladenosine (m6A) to mRNAs. While m6A-modified adenoviral RNAs have been previously detected, the location and function of this mark within the infectious cycle is unknown. Since the complex adenovirus transcriptome includes overlapping spliced units that would impede accurate m6A mapping using short-read sequencing, here we profile m6A within the adenovirus transcriptome using a combination of meRIP-seq and direct RNA long-read sequencing to yield both nucleotide and transcript-resolved m6A detection. Although both early and late viral transcripts contain m6A, depletion of m6A writer METTL3 specifically impacts viral late transcripts by reducing their splicing efficiency. These data showcase a new technique for m6A discovery within individual transcripts at nucleotide resolution, and highlight the role of m6A in regulating splicing of a viral pathogen.

Direct RNA sequencing on nanopore arrays redefines the transcriptional complexity of a viral pathogen.

Characterizing complex viral transcriptomes by conventional RNA sequencing approaches is complicated by high gene density, overlapping reading frames, and complex splicing patterns. Direct RNA sequencing (direct RNA-seq) using nanopore arrays offers an exciting alternative whereby individual polyadenylated RNAs are sequenced directly, without the recoding and amplification biases inherent to other sequencing methodologies. Here we use direct RNA-seq to profile the herpes simplex virus type 1 (HSV-1) transcriptome during productive infection of primary cells. We show how direct RNA-seq data can be used to define transcription initiation and RNA cleavage sites associated with all polyadenylated viral RNAs and demonstrate that low level read-through transcription produces a novel class of chimeric HSV-1 transcripts, including a functional mRNA encoding a fusion of the viral E3 ubiquitin ligase ICP0 and viral membrane glycoprotein L. Thus, direct RNA-seq offers a powerful method to characterize the changing transcriptional landscape of viruses with complex genomes.

RNA m6 A modification enzymes shape innate responses to DNA by regulating interferon ß.

Modification of mRNA by N6-adenosine methylation (m6A) on internal bases influences gene expression in eukaryotes. How the dynamic genome-wide landscape of m6A-modified mRNAs impacts virus infection and host immune responses remains poorly understood. Here, we show that type I interferon (IFN) production triggered by dsDNA or human cytomegalovirus (HCMV) is controlled by the cellular m6A methyltrasferase subunit METTL14 and ALKBH5 demethylase. While METTL14 depletion reduced virus reproduction and stimulated dsDNA- or HCMV-induced IFNB1 mRNA accumulation, ALKBH5 depletion had the opposite effect. Depleting METTL14 increased both nascent IFNB1 mRNA production and stability in response to dsDNA. In contrast, ALKBH5 depletion reduced nascent IFNB1 mRNA production without detectably influencing IFN1B mRNA decay. Genome-wide transcriptome profiling following ALKBH5 depletion identified differentially expressed genes regulating antiviral immune responses, while METTL14 depletion altered pathways impacting metabolic reprogramming, stress responses, and aging. Finally, we determined that IFNB1 mRNA was m6A-modified within both the coding sequence and the 3' untranslated region (UTR). This establishes that the host m6A modification machinery controls IFNß production triggered by HCMV or dsDNA. Moreover, it demonstrates that responses to nonmicrobial dsDNA in uninfected cells, which shape host immunity and contribute to autoimmune disease, are regulated by enzymes controlling m6A epitranscriptomic changes.

Molecular Aspects of Varicella-Zoster Virus Latency.

Primary varicella-zoster virus (VZV) infection causes varicella (chickenpox) and the establishment of a lifelong latent infection in ganglionic neurons. VZV reactivates in about one-third of infected individuals to cause herpes zoster, often accompanied by neurological complications. The restricted host range of VZV and, until recently, a lack of suitable in vitro models have seriously hampered molecular studies of VZV latency. Nevertheless, recent technological advances facilitated a series of exciting studies that resulted in the discovery of a VZV latency-associated transcript (VLT) and provide novel insights into our understanding of VZV latency and factors that may initiate reactivation. Deducing the function(s) of VLT and the molecular mechanisms involved should now be considered a priority to improve our understanding of factors that govern VZV latency and reactivation. In this review, we summarize the implications of recent discoveries in the VZV latency field from both a virus and host perspective and provide a roadmap for future studies.

Use of Whole-Genome Sequencing of Adenovirus in Immunocompromised Pediatric Patients to Identify Nosocomial Transmission and Mixed-Genotype Infection.

Background: Adenoviruses are significant pathogens for the immunocompromised, arising from primary infection or reinfection. Serotyping is insufficient to support nosocomial transmission investigations. We investigate whether whole-genome sequencing (WGS) provides clinically relevant information on transmission among patients in a pediatric tertiary hospital. Methods: We developed a target-enriched adenovirus WGS technique for clinical samples and retrospectively sequenced 107 adenovirus-positive residual diagnostic samples, including viremias (>5 × 104 copies/mL), from 37 patients collected January 2011-March 2016. Whole-genome sequencing was used to determine genotype and for phylogenetic analysis. Results: Adenovirus sequences were recovered from 105 of 107 samples. Full genome sequences were recovered from all 20 nonspecies C samples and from 36 of 85 species C viruses, with partial genome sequences recovered from the rest. Whole-genome phylogenetic analysis suggested linkage of 3 genotype A31 cases and uncovered an unsuspected epidemiological link to an A31 infection first detected on the same ward 4 years earlier. In 9 samples from 1 patient who died, we identified a mixed genotype adenovirus infection. Conclusions: Adenovirus WGS from clinical samples is possible and useful for genotyping and molecular epidemiology. Whole-genome sequencing identified likely nosocomial transmission with greater resolution than conventional genotyping and distinguished between adenovirus disease due to single or multiple genotypes.

Mode of virus rescue determines the acquisition of VHS mutations in VP22-negative herpes simplex virus 1.

It has been proposed that herpes simplex virus 1 with VP22 deleted requires secondary mutation of VHS for viability. Here we show that a replication-competent Δ22 virus constructed by homologous recombination maintains a wild-type (Wt) VHS gene and has no other gross mutations. By contrast, Δ22 viruses recovered from a bacterial artificial chromosome contain multiple amino acid changes within a conserved region of VHS. Hence, the mode of virus rescue influences the acquisition of secondary mutations.

Specific capture and whole-genome sequencing of viruses from clinical samples.

Whole genome sequencing of viruses directly from clinical samples is integral for understanding the genetics of host-virus interactions. Here, we report the use of sample sparing target enrichment (by hybridisation) for viral nucleic acid separation and deep-sequencing of herpesvirus genomes directly from a range of clinical samples including saliva, blood, virus vesicles, cerebrospinal fluid, and tumour cell lines. We demonstrate the effectiveness of the method by deep-sequencing 13 highly cell-associated human herpesvirus genomes and generating full length genome alignments at high read depth. Moreover, we show the specificity of the method enables the study of viral population structures and their diversity within a range of clinical samples types.

Leishmania-specific surface antigens show sub-genus sequence variation and immune recognition.

BACKGROUND: A family of hydrophilic acylated surface (HASP) proteins, containing extensive and variant amino acid repeats, is expressed at the plasma membrane in infective extracellular (metacyclic) and intracellular (amastigote) stages of Old World Leishmania species. While HASPs are antigenic in the host and can induce protective immune responses, the biological functions of these Leishmania-specific proteins remain unresolved. Previous genome analysis has suggested that parasites of the sub-genus Leishmania (Viannia) have lost HASP genes from their genomes. METHODS/PRINCIPAL FINDINGS: We have used molecular and cellular methods to analyse HASP expression in New World Leishmania mexicana complex species and show that, unlike in L. major, these proteins are expressed predominantly following differentiation into amastigotes within macrophages. Further genome analysis has revealed that the L. (Viannia) species, L. (V.) braziliensis, does express HASP-like proteins of low amino acid similarity but with similar biochemical characteristics, from genes present on a region of chromosome 23 that is syntenic with the HASP/SHERP locus in Old World Leishmania species and the L. (L.) mexicana complex. A related gene is also present in Leptomonas seymouri and this may represent the ancestral copy of these Leishmania-genus specific sequences. The L. braziliensis HASP-like proteins (named the orthologous (o) HASPs) are predominantly expressed on the plasma membrane in amastigotes and are recognised by immune sera taken from 4 out of 6 leishmaniasis patients tested in an endemic region of Brazil. Analysis of the repetitive domains of the oHASPs has shown considerable genetic variation in parasite isolates taken from the same patients, suggesting that antigenic change may play a role in immune recognition of this protein family. CONCLUSIONS/SIGNIFICANCE: These findings confirm that antigenic hydrophilic acylated proteins are expressed from genes in the same chromosomal region in species across the genus Leishmania. These proteins are surface-exposed on amastigotes (although L. (L.) major parasites also express HASPB on the metacyclic plasma membrane). The central repetitive domains of the HASPs are highly variant in their amino acid sequences, both within and between species, consistent with a role in immune recognition in the host.

COPASAAR--a database for proteomic analysis of single amino acid repeats.

BACKGROUND: Single amino acid repeats make up a significant proportion in all of the proteomes that have currently been determined. They have been shown to be functionally and medically significant, and are associated with cancers and neuro-degenerative diseases such as Huntington's Chorea, where a poly-glutamine repeat is responsible for causing the disease. The COPASAAR database is a new tool to facilitate the rapid analysis of single amino acid repeats at a proteome level. The database aims to simplify the comparison of repeat distributions between proteomes in order to provide a better understanding of their function and evolution. RESULTS: A comparative analysis of all proteomes in the database (currently 244) shows that single amino acid repeats account for about 12-14% of the proteome of any given species. They are more common in eukaryotes (14%) than in either archaea or bacteria (both 13%). Individual analyses of proteomes show that long single amino acid repeats (6+ residues) are much more common in the Eukaryotes and that longer repeats are usually made up of hydrophilic amino acids such as glutamine, glutamic acid, asparagine, aspartic acid and serine. CONCLUSION: COPASAAR is a useful tool for comparative proteomics that provides rapid access to amino acid repeat data that can be readily data-mined. The COPASAAR database can be queried at the kingdom, proteome or individual protein level. As the amount of available proteome data increases this will be increasingly important in order to automate proteome comparison. The insights gained from these studies will give a better insight into the evolution of protein sequence and function.